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Abstract 

A diagramatic heat kernel expansion technique is presented. The method is 
■ especially well suited to the small derivative expansion of the heat kernel, but 

it can also be used to reproduce the results obtained by the approach known as 
covariant perturbation theory. The new technique gives an expansion for the 
heat kernel at coincident points. It can also be used to obtain the derivative 
ffj ' of the heat kernel, and this is useful for evaluating the expectation values of 

the stress-energy tensor. 

Pacs numbers: 03.70. +k, 98.80.Cq 
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I. INTRODUCTION 



The heat kernel was first introduced into quantum field theory to help evaluate the 
polarization of the quantum vacuum caused by an external field |IJ . The heat kernel method 
invented by Schwinger was developed by DeWitt into one of the key techniques for evaluating 
the one loop effective action in the background field approach [@,[|. 

In many applications, the heat kernel is replaced by its asymptotic series in powers of 
the proper time |2|||,|[]. The coefficients in this series relate directly to physical quantities, 
with the first coefficient corresponding to the Weyl formula for the spectral density and 
further coefficients being relevant to the renormalisation group and to anomalies |J. The 
coefficients depend on external potentials, the field strengths of external gauge fields and the 
Riemannian curvature of the underlying space. The first six coefficients have been calculated 
by a wide variety of methods P,|7|- |10|1 - 

In practice, the proper-time expansion of the heat kernel is usually inadequate for ob- 
taining a good approximation to the effective action but an expansion in powers of covariant 
derivatives of the background fields may suffice. The terms in such a series can be organised 
in various ways depending on which quantities are to be regarded as being small. If the field 
strengths and derivatives of the external potentials are small then an expansion in powers 
of derivatives will be appropriate. This expansion is essentially a low energy expansion. 

Schwinger's original calculation of the vacuum polarization in a constant external field 
corresponds to a slightly different type of expansion where derivatives of the field strength 
can be neglected but the field strength itself may be large. This type of expansion leads 
directly to the Euler-Heisenberg effective action. 

A complimentary situation would be where the field strengths where small compared 
to their derivatives. The expansion should then be arranged in increasing powers of the 
curvatures with infinite numbers of derivatives resummed at each stage. This large derivative 
or high energy expansion is possible and it is used in an approach known as covariant 
perturbation theory [0,0 . 

In this paper we present a new diagramatic expansion technique, which is basically a 
small derivative expansion of the heat kernel, but which can also be used to reproduce the 
results obtained by covariant perturbation theory. The new technique gives an expansion 
for the heat kernel at coincident points. It can also be used to obtain local derivatives of 
the heat kernel, and this is useful for evaluating the expectation values of the stress-energy 



tensor [13[] 



There are close affinities between our method and the method invented by Barvinski 
and Vilkoviski and developed by Barvinski, Osborn and Gusev [[L4j]. In all of these 
methods, the diagrams are essentially the Feynman diagrams for a non-relativistic particle 
interacting with a background potential. However, there are important differences between 
the method presented here and previous methods for expanding the heat kernel. Firstly, 
we have organised our expansion into powers of derivatives of the potential, rather than the 
proper time. Our rules give the coincident heat kernel, wereas results in [O are given only 
for the integrated heat kernel. Furthermore, our rules apply in curved space, whereas results 
14f are for flat space, although a curved space generalisation was sketched out. 



in 



We will present various new results including new terms in the derivative expansion in 
flat space and for constant background gauge field strengths. We also give new terms in the 
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normal coordinate expansion of the metric that can be used elsewhere. 



II. THE BASIC METHOD 

Consider a Riemannian manifold (A4, g), and a second order operator A acting on a set 
of fields 4>{x). We are interested in the heat kernel G(x,x',t), defined by 

AG + G = 0, G(x,x',0) = 5{x,x') (1) 

where a dot denotes derivatives with respect to the proper time t and 5(x, x') is the covariant 
delta function. 

Throughout this section the operator will be taken to have a specific form 

- g^d^ v + + w (2) 

where <9 M is the partial derivative with respect to x^. The aim of this section is to obtain a 
series expansion for the heat kernel in powers of partial derivatives of g^, and W. 



Our method begins in the way introduced by Parker [I7| and follows some steps of 
reference ||15||. We define 



K(k,x',t) = J dfi{x) e~ ik{x - x,) G(x,x',t) (3) 
Taking the Fourier transform of equation (HD and integrating by parts gives, 

K(k, x', t)=-f dfi(x) {k 2 + W - Z) G{x, x', t) e - lk ^- x,) (4) 

Z(k,x,x') = (d„ - ik^){d v - ik v )h^ + (d„ - ik^ - W(x) + W(x') (5) 
and h^ v = g^ u — 5^ u . Expanding Z in a Taylor series about x = x' gives 

K(k, x', t) = -(k 2 + W - Z(k, D)) K(k, x', t) (6) 

where 

OO 1 

Z(k,D) = Y,-Z Mr D^...D^ (7) 

r=0 r ' 

and = id/dk^. 

Equation (^) has a solution in terms of a time-ordered exponential, 

K = e'^+^Texp (7* e^ 2+w ^'Z(k, D)e^ k2+W ^ dt^j (8) 

Redistributing some of the terms gives a more convenient form 

K = e- m e- fc2 '/ 2 Texp (j\ Wt 'Z{\i8',8)e- wt 'dS} e^ 2 (9) 



where 

8»(t') = -2ikH' + im (10) 

and 5' = —2ik. The function Z can be simplified by using the identity 

[F(6),5' ti ] = 2F(5) >fl , (11) 

which is valid for any function F(S), to remove some of the derivatives. After this simplifi- 
cation, 

Z=\hT{S)S'/ v + \T»{5)^ -W(6) (12) 
The coincidence limit of the heat kernel is obtained by integration 

G(x,x,t) = J d^(k)K(k,x,t) (13) 

where dfi(k) = d n k/(2n) n in n dimensions. Coincidence limits of derivatives of the heat 
kernel can be obtained in a similar way, for example 

[9fj,G] = -jdfi(k)^K(k,x,t) (14) 
[dfi'G] = J dv{k)K{k,x,t)% (15) 

where [. . .] denotes evaluation at x' = x. (The operators in equation @ act on the 5' in 
equation flT5p). Extra derivatives introduce extra factors of 5'. 
We shall define 

<->^/* We "' /2 - e "" 2 (16) 

where the normalisation factor 

K (t) = (AKty d ' 2 (17) 
for d dimensions. The coincident limit of the heat kernel is therefore 

G(x,x,t) =K (t)e- wt (Texp (J* e wt ' Z{\i8\8)e- Wt ' ' dA^. (18) 

Equation ( |TE| ) is the basic equation for obtaining the derivative expansion of the heat kernel. 
It yields a convenient diagramatic expansion by the conventional route using Wicks theorem. 
It is similar to an equation (for a more restriced operator) used in reference [|14 



We can expand the exponential and introduce creation and anihilation operators c±(k, D) 
that satisfy 

c^ e -fc 2 V2 = , [(£,<?+]= 5V (19) 
The operators S and 5' can then be replaced by c±, 
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n*o=y~{(*-*'K+*'4} (20) 

By Wick's theorem, the time-ordered product of operators can be replaced by products of 
propagators. There are three different propagators corresponding to the distinct ways of 
combining the operators, 

(T5»(U)5»(t J )) = 5^D(U,t J ) (21) 
(T# ll {t i mtj)) = VB(U,t j ) (22) 
(T6^t i )6' 1/ (t j )) = 6^D(t i ,t j ) (23) 

These can be evaluated using the creation and anihilation operators, 

D(t { , tj) = 2min(ti, tj) - 2r 1 t i t j (24) 
1$ (U, tj) = 29(t j - ti) - 2t-Hj (25) 
D (ti,tj) = -2r l (26) 

where 6{t) = 1 for t > and zero otherwise. 

After applying Wick's theorem, a typical term in the expansion of equation (|i~8|) will take 
the form 

K (t)e- Xt J dt 1 ...dt n T(e xtl Z 1 e- Xh . . . e xt " Z n e~ xt -) g^h . . .t n ) (27) 

where g(t± . . .t n ) is a combination of propagators and index contractions. The factors Zi 
arise from the Taylor series expansion of Z (see (|T2"D). When defining these factors it is 
convenient to multiply each term in the Taylor series by p\r\, where p is the order in 8' and 
r the order in 5, as this simplifies the combinatoric factors in the final expansion. 

The terms in the heat kernel expansion can be represented by diagrams, with lines 
representing propagators and vertices representing the terms Z^. There are three types 
of vertex, corresponding to the three terms in equation flT2|). Each term is divided by a 
numerical factor corresponding to the order of the symmetry group of the diagram in the 
usual way. 



III. CO VARIANT EXPANSIONS 

Second order operators are often written in the form 

A = -D 2 + X (28) 

where D = V + iA is a gauge covariant derivative, with gauge conection A and Levi-Civita 
connection V. The heat kernel expansion can then be expressed in terms of curvatures and 
covariant derivatives. We will use the notation 0(D n ) to describe the size of any term in 
the heat kernel expansion which has n covariant derivatives and any number of factors of 
X. The field strengths will be formally 0(D). 
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The partial derivative expansion can be obtained from the method just described. With 
this particular operator the vertices are generated by a function Z, where (see equation 

(0)), 

Z = \h^5'/ v + \T% - |z^{^, A v } - X 

-\i{hT, 5'^}A U - irM^ - sTApAv (29) 



We have made use of the identity (|il|) in this expression. The heat kernel expansion is 
generated by expanding equation fllSp. The anticommutators {, } only affect contractions 
at equal times, and this is taken into account in the diagrammatic rules given below. 

In principle, each of the terms generates a series of vertices when expanded in powers of 
5. However, most of the terms in Z can be eliminated if the propagator D attached to a 
vertex formed from the gauge field A is replaced by the propagator 

Dm (U, tj) = 26{U - t 3 ) - 2r 1 (30) 
The vertices are generated by 

Z = lh»%5l + \T% - Au } - X (31) 

This is particularly useful because the new propagator is the derivative of D . The details 
are given in appendix A. 

Each of the vertices in the diagramatic expansion corresponds to a multiple partial 
derivative of Z. It is convenient to define vertex functions by 

ikz e : ...ez.,,....,,,. (32) 

In a system of normal coordinates the vertex functions can also be expressed in terms 
curvatures and field strengths, as desribed in appendix B. 

The rules for evaluating the heat kernel G(x, x, t) consist of drawing all possible diagrams 
with lines from figure [l] and vertices from figure 0. The z'th vertex is associated with its own 
time variable ti, where < £j < t. The expression corresponding to the diagram has 

(1) For each line, stretching from a vertex % to vertex j, one of the propagators 

D(U, tj) = 2min(t;, tj) - 2r l Utj (33) 
~S(U, tj) = 29(t j - U) - 2r l tj (34) 
D(t i ,t j )=p8(t i -t j )-2r 1 (35) 

depending on the number of arrows. In these expressions, 6(0) = 1 and (3 = if only g and 
uj vertices are attached, whilst 9(0) = | and (3 = 2 if there is an A vertex. 

(2) For each vertex, the coefficient of the appropriate vertex function 

\e xti D\hre~ xti (36) 
i e xti D r T „ e -x ti ( 3? ) 

-ie xti D\A^- xti (38) 
-e xu D\Xe~ Xh (39) 
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These vertex functions can be expressed in terms of covariant tensors by using tables | to 



III. The components of these tensors are contracted according to the arrangement of the 



lines. 

The differentiated heat kernel at two coincident points is represented by diagrams with 
an external line for each derivative of the heat kernel. The free ends of the x and x' derivative 
lines are associated with the time variables t and respectively. (This respects the ordering 
of the terms in equation flT5|)). These contribute 

(3) For each external line, the appropriate propagator from rule (1) multiplied by | for 
an x' and — | for an x derivative. 

The resulting expression after applying rules (l)-(3) has to be time ordered and integrated 
with respect to the time variables. It is then multiplied by K exp(—Xt) (see equation fll7|)) 
and divided by the order of the symmetry group of the diagram. 

The expansion obtained from these diagrams is a derivative expansion where a diagram 
with n internal lines and m arrows produces a term of 0(D 2n ~ m ). 



IV. ABELIAN THEORIES IN FLAT SPACE 

The evaluation of the heat kernel simplifies considerably for abelian gauge theories in flat 
space. In some situations it is even possible to obtain exact results for the heat kernel, for 
example with a constant field strength F. In this section we will give the leading terms in a 
small derivative expansion of the heat kernel and show how resummations can be performed 
to generalise Schwinger's result for constant fields. We also show how resummations can be 
used to obtain an expansion for the heat kernel in powers of the curvatures, which can be 



regarded as a large derivative or high energy expansion []T2 . 

The most significant simplification for abelian theories is that only the connected dia- 
grams need be evaluated, since the results from the connected diagrams can be exponenti- 
ated to reproduce the results from the complete set of diagrams. This follows from the same 
simple counting argument that applies to the usual Feynman graph expansion. 

Other important simplifications include the fact that the vertex functions can be given 
explicitly, 

71 — 1 

D n f x A, = —-F v ^... Vn _^ . . . T"- 1 (40) 

/;r- v .v ; ,.,,r ...r (41) 

We can also reduce the number of diagrams by integration by parts. This is possible because 
of the relationships 

^(ti,t j ) = ^-D(t i ,t j ) (42) 

D(t i ,t j ) = -^~d(t i ,t j ) (43) 

that enable us to eliminate diagrams with D- 

The diagrams for some terms in the small derivative expansion with F = are shown in 
figure [| We write the small derivative expansion of the heat kernel in the form 
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G(x, x, t) = K (t) exp £W„ (44) 

\n=0 / 

Results up to 0(D 6 ) can be obtained by drawing all diagrams with up to three lines and 
using table [IV] to evaluate the time integrals, 

W = -Xt (45) 
W 2 = -l(d 2 X)t 2 + ±(dX) 2 t 3 (46) 
W A = -^(9 4 X)t 3 + ^(d^X^d^X)^ 

+ ±(d»Xd,d 2 X)t 4 - ±(d»X)(d»X)(d p d u X)t 5 (47) 
We = -^{d"X)t^ + ±.(d,d*XWX)t 5 

+^(d li d„d p X)(d»d»df>X)t 5 + ±(d y ,d v &X)(d»>9'X)1* 



(d„d 2 X)(d»d 2 X)t b - ^{d li X){d v d l X){dl i &'X)t b 



5040V ^ " 2520 

■±{d v X){d v X)(d»&'&X)t* - ^{d p X){d u d p X){d^d p X)t & 

■^{& l d v X){& , d p X){df , d^X)t % + ^(d^X)(d v X)(d p X)(d^d^X)t 7 

■ ggg (d,X) (d u X) (d^X) (d»d p X)t 7 (48) 



This exponentiated form also correctly reproduces any terms that can be reduced to products 
of 0(D e ) terms. As might be expected, most of the terms in W4 are also recovered in the 
proper time expansion of the heat kernel ||, or by other means [14|. However, most terms 
in Wq are new and were obtained with very little effort. 

The fact that the field strength has been regarded as a small quantity might be regarded 
as a drawback in the present approach. This can be overcome by a resummation of field 
strength terms to produce a small derivative expansion where F is of order D°. 

First of all define a new propagator Dp, 

00 

D F {t ri t s ) = Y,D n {t r ,t s )F n (49) 

n=0 

where F is a matrix and 

D n (t r , t.) = j dt x ... dt n D(t r , ti)^(ti, t 2 ) • • • Q (50) 

We also define Dp by differentiating with respect to t r . An explicit expression for Dp and 
is given in appendix C. 

The contribution to the heat kernel from ring diagrams involving only the field strength 

is 

Wo = E 7T tT ( Fn ) / ^n-l(t',t')dt' (51) 
2n Jo 

We can write this in terms of the new propagator, 

W = \ f F du f dt'^^t'^t') (52) 
2, Jo Jo 

By equation ( p5|) , 
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1 i sinhFt , . 

W = --trbg— — (53) 

This recovers Schwinger's result for the heat kernel in a background with constant field 
strength fl]. 

In other diagrams, replacing the propagator D by Dp resums all of the terms involving 
F. For example, 

W 2 = \J dt x dt 2 X^X^D F {t x ,t 2 f v - \ J dt 1 tiX )lw D F (t l ,t 1 )' a/ . (54) 
Again, after using equation (|U5|), 

W 2 = (\X tlx X >v - \X tlw ) (F-H 2 cothFt - F-Hy (55) 

This new result reduces to the previous result if F — 0. 

A quite different situation exists when the derivatives are large, d 2 X ^> X 2 . The large 
derivative expansion of the heat kernel is then given by the diagrams shown in figure f|. By 
summing these diagrams it is possible to order the derivative expansion in powers of X, thus 

G(x, x, t) ~ Gx(x, x, t) + G2(x, x, t) + . . . (56) 

where 

Gi(x,x,t) = K (t) J dfiik^i-^e^fiik^tXih) (57) 

G 2 (x,x,t) = K (t) J c/ / u(A; 1 )rf/i(A; 2 )e* {fcl+fc2) - a: /(^it,^ 2 t)t 2 X(A; 1 )X(A; 2 ) (58) 

We have written X(k) for the Fourrier transform of X. According the the diagrammatic 
rules, the single vertex diagrams give 

Wit) =£ll| *lD(t u hYiih) 2 ? (59) 

Similarly, for the two-vertex diagrams 

f(ik 1 t,ik 2 t)= (60) 

E J ^^^(*i,«i)^(*i,«aW* a ^) r (i*i) a, '(i*2) ap (-*i " 

Using the explicit form of the propagator (|33D , we can perform the sumation to express f(a) 
and f(a, b) in closed form, 

f(a) = f 1 dxe x(1 ~ x)a2 (61) 
Jo 

f(a, b) = 2 [ dxx [ Xl dx 2 e Xlil -' Xl)a2+X2{1 ~ X2)b2+2xiX2a - b . (62) 
Jo Jo 



The local expression for the heat kernel in the large derivative limit given by equations 
and can be used in any region in which the derivatives are large, even if they are not 
large everywhere. 

For the integrated heat kernel, k\ = —k 2 and we can use f(a, —a) = /(a) to obtain 



dii{x)G{x,x,t) = K (t) I dfi(k)f(ikt)t 2 X(k)X(-k) (63) 
This demonstrates agreement between the average of our heat kernel and the results for the 



integrated heat kernel given in reference [JL2 
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V. CONCLUSIONS 



We have seen that a diagrammatic approach gives the covariant derivative expansion 
of the heat kernel in both the large and small derivative limits. The method gives the 
coincident limits of the heat kernel or derivatives of the heat kernel and can be applied to 
any minimal second order operator. 

The new diagramatic expansion reproduces known results quickly and reliably, and we 
have also obtained some new terms in the derivative expansion in flat space with constant or 
vanishing background gauge field strengths. The derivation of more new results is presently 
underway. 

One important area of application involves finding the heat kernel at finite temperatures 
to evaluate thermodynamic functions. The new method generalises very easily to this situ- 
ation by following the ideas presented in ||15|| . If only spatial derivatives are present, then 
the only change is the replacement of K Q (t) by a Jacobi theta function, 

K (t) = (47rt) (1 ^ )/2 # 3 (0,47u/T 2 t) (64) 

where j3 is the inverse temperature. 



APPENDIX A: VERTEX ELIMINATION 

In this appendix we examine how some of the vertices in the diagramatic expansion can 
be eliminated by modifying the propagator. 

The full set of terms which generate the vertices are given by equation (f29|), 

z g = Ihrs^i, z^ = |r% z a = -\i8^%, A u }, 

Z ag = -\ih^{5'^A u }, Z ua = -ir»An, (Al) 
Z a a = ~ A^Ay, Z aga = —h^A^Ay. 

Modified brackets (. . .) m are defined by applying Wick's theorem with a modified propagator 

D= 25(U - tj) - 2t~ l (A2) 

connecting the A vertices and with a restricted set of vertices generated by Z g , Z w and Z a 
only. Contractions between the 8' factors reproduce the missing vertices. Expanding the 
exponential in equation (^j) with the reduced set of vertices gives a typical term 

1 



p\q\r 



■(TZ*Z«Z:) m = £ c iUm (TZi g Zt a Z l M Z™ a ZlZlZi) (A3) 



where 

p' = p — j — m, (A4) 

q' = q - k, (A5) 

r' = r — j — k — 21 — 2m (A6) 

and the Cjkim are combinatorial factors. These factors can be found by counting the number 
of ways in which to divide up the vertices and then counting the possible contractions. The 
result is equal to the coefficient of the identical term obtained by expanding the time ordered 
exponential with the full set of vertices and the unmodified propagator. 
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APPENDIX B: NORMAL COORDINATE EXPANSIONS 



It is well known that in a system of normal coordinates the partial derivatives of the 
metric and the gauge fields can be replaced by covariant expressions involving curvatures 



17fl . We will give a brief review here to bring out some important features and to give the 
results in a form that can be used in the diagramatic expansion of the heat kernel. 

The gauge covariant derivative along a coordinate basis will be denoted by D^. The 
curvature operator 1Z can be defined by the covariant derivative acting on a vector field e a , 

[£>„, D v ]e a = K(e^ e u )e a (Bl) 

By allowing the vector e a to point in either the direction tangential to the manifold or the 
direction of the internal symmetry allows the connection coefficients to include both the 
tetrad connection and the gauge field A^. The corresponding components of the curvature 
operator are the Riemann tensor and the Field strength tensor. 

We choose a point P and an orthonormal frame e a at P to set up the normal coordinates. 
Consider a family of geodesies passing through P with tangent vectors £. The normal 
coordinates x M of a point Q are defined by 

x » = o^(t) (B2) 

where a is the distance along the geodesic from P to Q and £ (r) is the unit tangent vector at 
P. This construction implies that the coordinate vectors e M satisfy a commutation relation 

Ds(ae„) - oD^ = A£ (B3) 

for a constant A. We can also arrange that D^e a = for the reference frame used to calculate 
the curvature components. 

By repeatedly differentiating aA^ and ae^ at P we get two important identities 

71 — 1 

£>r A » = —DT'nt, e M ) (B4) 

71 — 1 

Dte ll = ——D£- 3 R(Z,e ll )Z (B5) 

5 71+1 s 

These can be used recursively to allow us to replace the connection components and their 
derivatives with covariant derivatives of the field strengths. 

If the manifold is flat then equation (|B4|) has an explicit solution 



77 — 1 

D% Ap = F m . Mn ^r • • • r- 1 (B6) 

where the ';' denotes the gauge covariant derivative. In this case the choice of gauge is 



identical to the gauge introduced originally by Fock and developed by Schwinger [pL8| , pL9 . 

In the curved case, we are interested in derivatives of the inverse metric and the connec- 
tion T M , where 



= g» u uj >v - W\ v (B7) 
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and uj = log \g^ v \ l l 2 - We can use matrix notation, writing e for the tetrad and g for the met- 
ric, with g = e T e and u = — log |e|. Covariant expressions for derivatives of the determinant 



and the inverse of the metric at P can then be obtained from Faa di Bruno's formula |20 
for example 

n (-l) m r?l 

D n u =T T - —7-T, ~7 — 1\ — {tr ((De) ai . . . (D n e) a ") + permutations} (B8) 

D n g- X = E £(-l) m a ,L 7 |w • • • (IPe)** + permutations} (B9) 

summed over a± + a 2 + . . . a n = m and a\ + 2a 2 + . . . na n = n. Expressions for D n e can be 
obtained from equation (|B5| ). 

Our results are tabulated in tables | to TT|. These can be partially checked against 
similar expansions in a paper by van de Ven ||, which contains results up to n = 6. Our 
Rj corresponds to |±|Kj +2 and F,- to £jj|Yj- + i in his 'index-free' notation. 



APPENDIX C: TIME INTEGRALS 



Each diagram T> in the heat kernel expansion in derivatives corresponds to an expression 
of the form 



G v = K (t)e- Xt J dt x ... dt n T (e xtx Z x e~ x ^ . . . e xtn Z n e~ xtn ) gfa . . .t r , 



(CI) 



We will consider the abelian case first. 
When the terms commute we have 



G v = K {t)e~ xt Z l ...Z n Jdh... dt n g{t x ...t r 



(C2) 



where g(t± . . . t n ) is a product of propagators. The integrals can be reduced recursively by 
equations such as 



* D(t 1 ,h) p dt 2 = -!—D(ti,tiY 
o p + l 



(C3) 



or similar results. We can also replace arrows over the propagators D(ti,tj) by derivatives 
with respect to U for — > and tj for <— . These simplifications reduce the integrals for terms 
up to 0(-D 6 ) to those tabulated in table [TV . 

The results quoted in section 3 can be obtained from the fact that D w (ti, tj) satisfies 



d 2 D u 
dt) 



2u 



dD u 



25^ -U) 



(C4) 



and D w vanishes at and t. Consequently 



D w (t{ , tj / 



-l)(e- 



-2vt\ 



( e 2w(t,-t) _ l)( e -2urt, _ X ) 
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ti ^> tj 



to *\ tj 



(C5) 



Also D u (U,tj) = diD u (ti,tj). 

In the non-abelian case can need to introduce a sum over permuations ir, 

G v = K (t)e- Xt F v (Lit)Z v{1) . . . Z <n) (C6) 
where Lj is the commutator LjZj = [X, Z^Sy and 

F n ( ai ) = I dh... dt n g(h . . . tn )e^ + - a ^ t (C7) 

■'*7r(l)>— *ir(n) 

The functions F n depend on the diagram but not on the particular gauge group. 

We can process the results further by using a gauge transformation to diagonalise the 
matrix X. It is then possible to choose a canonical basis E a of the Lie algebra with the 
property that 

[X,E a ] =w(a)E a (C8) 

for constants w(a). The vertices which do not commute with X belong to representations 
of the Lie algebra of the gauge group, Zj = Z°;E a , therefore 

G v = K (t)e- Xt ^ c ai ... an Z? . . . Z? (C9) 

ai...a n 

where 

c ai ... an =J2F 7r (w(a i )t)E an{1) . . . E a<n) . (CIO) 

7T 

We see from this expression that the range of values where F n needs to be calculated is 
restricted by the Lie algebra roots. 
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TABLES 



D n g~ 1 



2 — |i?o 

3 -Ri 

4 -p2 + §*§ 

5 — 3^3 &RqRi 

c 10 n i 100 p p I 85 p2 160 p3 

7 3d 1 122 p p 1 357 p p 1 6611 p2 p 1 661 p p p 

I — 2-K5 + —^(0-^3) + —■ K (l- K 2) + "IfT^O^l) + -g--Ko-Kl-Ko 
o 14 p 1 280 p p 1 824 p p 1 308 p 2 1596 p2 p 

8 —-q-Kq + -g-Vt(o-«4) + —^(1-^3) + ~ K 2 — R( R2) 

3052 r>, p2 868 p p p 4676 p p p , 628 p4 
g - -K(0- K i) g--f£i/£oili g— ito-ri2-rto + — -flQ 



TABLE I. Expressions for 1 in terms of curvatures. In this table, Rj denotes the linear 

mapping RjX = (DlTl)(£,X)£, where 1Z is the curvature operator. In terms of components, 



D n g »v = e ^ . ( D n g -r )ev and e ^ . R 3 e u = R^ lvm .. Uj+2 ^ . . . 



n expression 

2 — |i?o 

3 — T^Ri 

4 -5-R2 + ^-Ro 2 

5 — 3-R3 + ^RqR\ 

6 5 p i8pp 1 15 p 2 16 p 3 
— yit4 + + 14-n-l — 63^° 

7 — § -R5 + 4:RqRs + 5-R1-R2 — -^Ro^i 

o 7 p 1 20 p p 1 70 p p 1 28 p 2 272 p 2 p 310 p p 2 , 16 p4 

° — 9-K6 + -9--K0-K4 + -9--K1-K3 + "5"-K2 — ^5"-K0 -^2 9--K0-K1 + JE K 



TABLE II. Expressions for D^u, where uj = log \g tJlV \ 1 ^ 2 , are obtained by tracing the entries in 
this table. 



D n - l A 



2 \Fq 

3 pi 

4 4-F2 + 3-F0-R0 

5 I-F3 + I-F1-R0 + \FqR\ 

6 + l-F2.R0 + \F\Ri + \FqR 2 + |Fbi?o 

7 yi^ + 'y-FzRo + ^^i?! + ^-F\R2 + jFiRq + 7.F0.R3 + ^FqRiRo + ^FqRqRi 

S 'p±^pp i35ijp i63pp i21r p2 1 7 p p 1 7 77 p p 

o g^6 + "g-^4-Ko + X^ 3 ^ 1 + T 2_K2 + T 2jK + ^-^i-^ 3 + ^-^l^l^O 



+ 4-Fii?o^?l + §Fbi?4 + ^-FqR 2 Rq + ^-FqR\ + ^FqRqR 2 + ^0^0 

TABLE III. Expressions for _D^ -1 v4 in terms of curvatures and field strengths. In this table 
Fj denotes the linear mapping FjX = (D^TZ)(^,X) with gauge field strength indices. In terms of 
components, D^A^ = (D^A)e^ and F^ = F vm .„ Vj+1 ^ . 
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vertices 


9 


numerical factor 


term 


1 


D{hM) p 


2P( P \) 2 /{2p + l)\ 


t P+i 




D(t 1 ,t 1 )P^(t 1 ,t 1 ) 2 i 


2Pp\{2q)\{p + q)\/(2p + 2q + l)\q\ 


t P+l 


2 


D(h, h)D{t u t 2 )D(t 2 ,t 2 ) 


17/630 


t 5 


n 


D(t 1 ,t 2 )...D(t n ,t 1 ) 


{-2f n - l B 2n /{2n)\ 


f 2n 


2n 


'fi{t 1 ,t 2 )...r${t 2n ,t 1 ) 


-2 2n - 1 B 2n /{2n)\ 


f 2n 



TABLE IV. The time integrals corresponding to diagrams for an abelian operator up to 0(D 6 ) 
can be reduced to the examples tabulated here. The results are given by the numerical factors 
multiplied by the terms in t, and B n are Bernouilly numbers. 
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FIGURES 



D(tiJj) D(t iy tj) D(ti.tj) 

FIG. 1. Lines appearing in the expansion of the heat equation. They correspond to the 
two-point functions of time as indicated. 



\ r \ r \ 

g) f (»)— ► f ( A ) ► F (X 







\D r gLiv \Dr<i)V -iD^A/j, DrX 

FIG. 2. Vertices appearing in the diagrammatic expansion. These vertices are attached to r 
plain lines and up to 2 lines with arrows. Each vertex corresponds to an integral over time and the 
vertex function indicated below the vertex. Covariant expressions can be found in the tables. 

Q?D (Tt=© ©— (*D 



FIG. 3. The terms in the derivative expansion of 0(D 4 ) are formed from the diagrams shown 
here. 



x V fix] = 





FIG. 4. Diagrams of this type dominate the large derivative approximation. 
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